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' Abstract. We study the resonant absorption of MHD waves in magnetized flux 

' tubes with a radial density inhomogeneity. Within the approximation that resistive 

and viscous processes are operative in thin layers surrounding the singularities of 



Or 



O , 

analyzed 



the MHD equations, we give the full spectrum of the eigenfrequencies and damping 
rates of the MHD quasi modes of the tube. Both surface and body modes are 
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1. Introduction 

Ionson (1978) was first to suggest that the resonant absorption of MHD waves in coronal 
plasmas could be a primary mechanism in coronal heating. Since then, much analytical 
and numerical work has been done on the subject. Rae and Roberts (1982) investigated 
both eikonal and differential equation approaches for the propagation of MHD waves 
in inhomogeneous plasmas. Hollweg (1987a, b) considered a dissipative layer of planar 
geometry to study the resonant absorption of coronal loops. Poedts et al. (1989, 1990) 
developed a finite element code to elaborate on the resonant absorption of Alfven waves 
in circular cylinders. 

Davila (1987) and Steinolfson & Davila (1993) did much work on resonant absorption 
through resistivity. Ofman et al. (1994) included viscous dissipations in their analysis 
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and concluded that the heating rate due to shear viscosity is comparable in magnitude 
to the resistive resonant heating. Also, they concluded that the heating caused by com- 
pressive viscosity is negligible. Goossens et al. (2002) used the TRACE data of Ofman & 
Aschwanden (2002) to infer the width of the inhomogeneous layer for 11 coronal loops. 
Rudcrman & Roberts (2002) did similar analysis with the data of Nakariakov et al. (1999). 
Van Doorsselacrc et al. (2004a) used the LEDA code to study the resistive absorption 
of the kink modes of cylindrical models. They concluded that, when the width of the 
nonuniform layer was increased, their numerical results differed by as much as 25% from 
those obtained with the analytical approximation. Van Doorsselaere et al. (2004b) inves- 
tigated the effect of longitudinal curvature on quasi modes of a typical coronal loop. They 
found that the frequencies and damping rates of ideal quasi modes were not influenced 
much by the curvature. Andries et al. (2005) studied the effect of density stratification 
on coronal loop oscillations, and conclude that longitudinal mode numbers are coupled 
due to the density stratification. 

In the absence of resonance, Edwin & Roberts (1983) and Roberts et al. (1984) 
gave a comprehensive account of the theoretical and physical concepts related to coronal 
oscillations. Karami et al. (2002; hereafter paper I) studied the full spectrum of MHD 
modes of oscillations in zero-/? magnetic flux tubes with discontinuous Alfven speeds at 
the tube's surface. In the vicinity of singularity, field gradients are large. Recognizing 
this, Sakurai et al. (1991a, b) and Goossens et al. (1992, 1995) developed a method to 
analyze dissipative processes in such regimes and to neglect them elsewhere. 

Here we combine the two techniques of paper I and of Sakurai et al. (1991a) to obtain 
the resonant damping rates for the full spectrum of the normal modes of magnetic flux 
tubes. 



2. Equations of motions 

The linearized MHD equations for a zero-/3, but resistive and viscous, plasma are 

^ = -!-{(V x SB) x B + (V x B) x SB} + ^V 2 5v, (1) 
at Aitp p 

? = Vx(kxB) + V 2 (5B, (2) 

Ot 47TCT 

where Sv and <5B are the Eulerian perturbations in the velocity and the magnetic fields; 
p, a, r], and c are the mass density, the electrical conductivity, the viscosity and the speed 
of light, respectively. The simplifying assumptions are: 

— Under coronal conditions gas pressure is negligible (zero- (3). 

— Density scale heights are much larger than the dimensions of flux tube, so that the 
gravity stratification is negligible. 

— The tube geometry is a circular cylinder with coordinates (r,<fi,z). 
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— There is a constant magnetic field along the z axis, B = Bz. 

— The equilibrium is static. 

— There is no initial steady flow inside or outside of the tube. 

— Viscous and resistive coefficients, 77 and a respectively, are constants. 

For a variable density, p(r), a singularity develops wherever the local Alfven frequency 
becomes equal to the global frequency of the mode. The relevant radial wave number 
vanishes and resonant absorption takes place. Let us denote the radius of the tube by 
R and a radius beyond which the resonance occurs by R\ < R. The thickness of the 
inhomogeneous layer, a = R — R\, will be assumed to be small and will be arbitrarily 
taken to be of the order R/10. The choice of density profile is also unimportant. We will 
assume two constant densities, pi in r < R\ and p e < pi in r > R, interconnected with 
a linearly varying profile in R\ < r < R. 

In the remainder of this section the following steps are taken: 

a) In r < Ri and r > R, dissipative terms are neglected. Solutions of Eqs. JU and 
are obtained as per paper I, and their differences across the inhomogeneous layer are 
calculated. 

b) In i?i < r < R, within which the resonant layer resides, solutions are found by 
expanding Eqs. Q and J5J) around the singular point, and the jumps across the 
resonant layer are found by the prescriptions of Sakurai et al. (1991a). 

c) Interior and exterior solutions are connected by requiring the jumps in (b) to be equal 
to the differences in (a). This gives an analytical expression for a complex dispersion 
relation to be solved for the frequencies and the damping rates. 



2.1. Interior and exterior solutions 

From paper I, in the absence of dissipations, all components of Sv and the transverse 
components of are expressible in terms of SB Z only. The latter, in turn, is the solution 
of a second order differential equation. Thus, 



r dr 
5B r = 



dSB z 



k 2 dr 
ik z dSB z 
k 2 dr 



+ (k 2 - — )5B Z = 0, 



SB, 



mk 7 SB, 



(3) 
(4) 



i uj d8B z m uj 8B Z 

SVr = -k^B^r-' = PflV' * W « = « (5) 



where k 2 — uj 2 /v\ — k 2 , and VA(r) — B j \J Anp(r) is the local Alfven speed. Here, we 
have assumed an exponential </>, z, and t dependence, exp[i{m<j) + k z z — ojt)] for any 
component of Sv and (5B. 
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In the interior region, r < R\, solutions of Eq. © are 

Im(\ki\r), ki < 0, surface waves, 

8B Z — ^ J m (\ki\r), kf > 0, body waves, (6) 

kf =ui 2 /v\ z -k 2 z , 

where J m and I m are Bessel and modified Bessel functions of the first kind, respec- 
tively. In the exterior region, r > R, the waves should be evanescent. Solutions are 

8B Z = K m {k e r), k 2 e = k 2 z - u 2 /v 2 Af > 0, (7) 

where K m is the modified Bessel function of the second kind. 

2.2. Dispersion relation and damping 

From Sakurai et al. (1991a), Goossens et al. (1992, 1995), and Erdelyi et al. (1995), the 
jump conditions across the boundary for 5B Z and 5v r are 



[Sv r ] =~™— 2 B z SB Zl (9) 



[SB,] = 0, (8) 
1 

|A|pMrV 

where R\ < ta < R is the radius at which the singularity occurs and k 2 {rA) = 0, 
u = w + 17, and A = -B 2 ^{y)\r A - Substituting the fields of Eqs. ©, © and (J7J 
in jump conditions and eliminating the arbitrary amplitudes of the waves, as foreseen 
initially inside and outside of the boundary layer, gives the dispersion relation 

do(uj, m) + di(ui, m) = 0, (10) 

where 

_ l K' m (\k e \R) , l J'„A\h\Ri) hU 

a ^^- ke Km (\k e \R) + h JMRiV [ ' 



di{Cj,m) = -in— — ^ 2 . (12) 



1 

\&\p(rA)r 2 A 

In principle tl> = uj + i"/ is expected to be found as a solution of Eqs. l|lUll2f> . In particular 
for 7 <c oj, Eq. I|1U|) can be expanded about ui to give 
Im(di(u), m)) 



(13) 



dd (u,m)/dQ - =w ' 

The results for surface waves are the same as Eqs. (jlOII - <|13[) . except that J m is replaced 
by I m everywhere. This completes the formal solutions of ui and 7. Further analytical 
progress is still possible when the inhomogeneous layer is thin enough. 



3. Thin boundary approximation 

Here we assume (R — R\)/R = a< 1. Equation ijTUIl reduces to d ~ 0. For body waves, 
the latter is studied in ample detail in paper I. There, a trio of wave numbers (n, m, I), 
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corresponding to r, (j> and z directions, is assigned to each mode. From Eq. (|13[) the 
corresponding damping rate becomes 

d f 1 K' m (\k e \R) 



lnml i ^l ml (p l -p e )R^j / ~dw \ k e K rn {\k e \R) 

i J' m {\h\R) \ (u] 
hJ m {\h\R)]- 1 ' 

Again the results for surface waves are the same as those for body waves except that J m 
is replaced by I m . The two waves exhibit differences, for J m and I m behave differently 
at the boundary. We also note that each surface mode is designated by only two wave 
numbers, (m, I) corresponding to (f> and z directions, respectively. 

For m = 0, resonant absorption does not take place, because the jumps in 5B Z and 
Sv r vanish and Eqs. 110H13f) are not valid anymore. 



4. Numerical results 

As typical parameters for a coronal loop, we adopt radius = 10 3 km, length = 10 5 km, 
Pi = 2 x 10~ 14 gr cm~ 3 , p e /pi — 0.1, B = 100 G. For these parameters one finds 
— 2000 km s _1 , va c — 6400 km s -1 and u>a = 2 rad s _1 . 

In Fig.nthe magnetic field component, 8B z (r) 1 for n = 1, m = 1 and I = 41, 50, 100, 
are plotted versus x = r/R. They are normalized to max(SB z ,l = 41). The segment of 
the plot in < x < 1 is from the interior solutions of Eq. JfjJ, while the segment in x > 1 
is from the exterior solutions of Eq. Q. Slopes are discontinuous at x = 1. Actually, the 
correct location of the discontinuities is the point of singularity, ta- Its appearance at 
x = 1 is an artifact caused by extrapolating the interior solutions to cover the boundary 
layer, rather than using the exact solutions there. As I increases, the maximum wave 
amplitude moves towards the tube axis and away from the inhomogeneous layer. 

In Fig. [3 the frequencies and the damping rates are plotted versus I. The frequencies 
increase with increasing n and/or /. For n = 1 and 2, damping rates exhibit maxima at 
/ w 50 and 70, respectively. For n = 3 there is only a declining branch towards higher I 
values. Our explanation for the occurrence of maxima in damping rates is the localization 
of the maximum amplitude of a wave within the resonant layer, where the dissipation is 
expected to be the highest. As Fig. ^ shows, this maximum moves away from that layer 
at both lower and higher / values. 



For surface waves, specified by two mode numbers (m,l), u m i and 7,„/ are plotted in 
Figs. 01 and Both frequencies and damping rates show a monotonous increase with I. 
Two tube radii, R/L = 0.01, & 0.02, are considered here. The frequencies and damping 
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Fig. 1. Magnetic component 5B z /5B zmax ^ i=41 , for body modes, where n = 1, m = 1 
(kink), and I = 41,50,100. Slopes are discontinuous at x = 1. The maximum wave 
amplitude shifts towards the tube axis as I increases. 
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Fig. 2. Body modes: u} nm i (solid lines) and 7 nm ; (dashed lines) versus I for m = 1 
(kink) modes. Auxiliary parameters are: the tube radius 10 3 km, the tube length 10 5 
km, B = 100 G, pi — 2 x 10~ 14 grcrn -3 , p e /pi = 0.1. Frequencies are in units of the 
interior Alfven frequency, uja — 2 rad s . 




rates of the thicker tube are almost double those of the thinner tube. This is expected, 
because in the limit of thin tubes, both oj and 7 are proportional to R/L (see, e. g., Van 
Doorsselaere et al. 2004) . The frequencies of the surface waves of Fig. and the n = 1 
body frequencies of Fig.|21both behave similarly. This is because of the similar behavior of 
Im{kr) and J m (kr) at small arguments. In physical terms this means, at least in the thin 
tube approximation, the surface and the r- fundamental body waves contribute in equal 
manner to the heating mechanism; they both dump the bulk of their energies near the 
tube surface, and perhaps both are excited with comparable amplitudes and energies. The 
w/27T7 is the number of oscillations taking place before a wave is completely attenuated. 
For thick and thin tubes, and 1 < I < 10, this number is about 1.7 and 8, respectively. 
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An observed value of Nakariakov et al. (1999) from TRACE data is about 3.39. Those 
of Goossens et al. (2002) from the same source range from 1.01- 3.21 corresponding, 
according to the authors, to various values of R/L, a/R, and p e / Pi- 
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Fig. 3. Surface modes: w m i and 7 m /versus I for rn = 1, R/L = 0.01 (Solid line) & 0.02 
(Dashed line). Auxiliary parameters as in Fig. [3 
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Fig. 4. Surface modes: u> m i and 7 m / for m = 2 and R/L = 0.01. Auxiliary parameters as 
in Fig. El 



5. Summary 

We studied the MHD quasi modes of coronal loops. On the assumption that ohmic and 
viscous dissipations are operative within a thin boundary layer, we obtained an analytic 
dispersion relation and solved it numerically for the mode frequencies and the damping 
rates. For realistic values of the initial parameters thickness- to- height ratio of the loop, 
the density contrast with the background medium, and the equilibrium magnetic field- 
our numerical values agree with those obtained from observations. As the longitudinal 
wave number increases, the maximum amplitude of the body eigcnmodcs shifts away from 
the resonant layer and causes a decrease in damping rates. Its behavior with increasing 
radial wave number is not, however, all that straightforward. In additional, we have 
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shown that body and surface modes may contribute equally to the heating of coronal 
loops. 
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